Overlimiting current in a microchannel 
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We revisit the classical problem of diffusion-limited ion transport to a membrane (or electrode) by 
considering the effects of charged side walls. Using simple mathematical models and numerical sim- 
ulations, we identify three basic mechanisms for over-limiting current in a microchannel: (i) surface 
conduction carried by excess counterions, which dominates for very thin channels, (ii) convection 
by electro -osmotic flow on the side walls, which dominates for thicker channels and transitions 
to (iii) electro-osmotic instability on the membrane end in very thick channels. These intriguing 
electrokinetic phenomena may find applications in biological separations, water desalination, and 
electrochemical energy storage. 

PACS numbers: 



Introduction. — All electrochemical cells involve ion- 
selective surfaces, such as electrodes or membranes, that 
pass current carried by certain active ionic species [T]. 
In an unsupported bulk electrolyte, the rejection of the 
opposite species (to maintain neutrality) leads to salt de- 
pletion and "concentration polarization" (CP). Theoret- 
ically, at the diffusion- limited current, the solution con- 
ductivity vanishes at the surface, but in practice, "over- 
limiting current" (OLC) is often observed and associated 
with an extended depletion zone (Fig. la). 

Possible mechanisms for OLC have been debated for 
decades. In the context of water electrodialysis [2J, it 
is now understood that OLC can be due to chemical ef- 
fects (water splitting or exaltation), which create addi- 
tional ions, or a physical effect, electro-osmotic instability 
(EOI) 3 J , which enhances ion transport by convection [3] 
(Fig. Id). The flow is driven by extended space charge 
(ESC) in the non-equilibrium double layer on the selec- 
tive surface [5|. Similar phenomena can also occur in 
dendritic electrodeposition [7j , electrochromatography [5] 
and capacitive desalination at large voltages [5]. 

Some experiments have shown that CP is also af- 
fected by geometrical confinement. In thin-gap (100/im) 
electrodeposition cells, gravitational and electro-osmotic 
convection can affect the spreading depletion layer [7]. 
In thinner microfluidic devices, steady depletion layers 
in transverse flows can form surprisingly sharp inter- 
faces ("demixing"), whether triggered by capacitive de- 
salination [9] or electrodialysis [10] . Additionally, convec- 
tion can be suppressed by decreasing the channel height, 
thereby reducing (but not eliminating) OLC [TTJ [T2J . In 
very thin (1/xm) channels, it was recently shown that CP 
interfaces can propagate at constant current, analogous 
to shock waves in gases [T3] ■ These "desalination shocks" 
result from geometry-dependent surface conduction (Fig. 
lb) [14], but the role of convection and the connection 
with OLC have not been considered. 

In this Letter, we develop a unified theory of steady CP 
in a microchannel (neglecting chemical effects) and iden- 
tify physical mechanisms for OLC related to charged side 




FIG. 1: Physical picture of over-limiting current (OLC) in 
a microchannel from a reservoir (left) to a cation-exchange 
membrane (right). The volume-averaged conductivity profile 
(a) exhibits a classical linear diffusion layer (continuous line) 
and a constant depleted region (dashed line), where the cur- 
rent is carried primarily by (b) SC, (c) EOF or (d) EOI with 
increasing channel thickness. 



walls. We show that surface conduction (SC) dominates 
in very narrow channels (Fig. lb), while electro-osmotic 
flow (EOF) on the side walls dominates in wider channels 
(Fig lc) and transitions to EOI on the selective surface 
in very wide channels (Fig. Id). We are not aware of any 
prior mention of the SC mechanism for OLC, whereas our 
theory of the EOF mechanism differs from the recently 
proposed ID Taylor- Aris dispersion model (strictly valid 
for aspect ratios > 10 3 and zeta potentials > 1 V) [16], by 
taking into account the 2D velocity field at the channel 
end. 

Our analysis is based on the following model problem. 
A symmetric, binary z : z electrolyte fills a parallel-plate 
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(2D) microchannel of length L and thickness H, which 
is open to a well-mixed reservoir of concentration cq at 
x = 0. A voltage V is applied across the microchannel to 
drive a current I through an ideal cation-selective solid 
surface (the "membrane") at x — L, which could repre- 
sent a cation-exchange membrane, a negatively charged 
nanochannel, or a metal electrode undergoing electro- 
plating. The side walls have a charge density <j s < 0, 
which we take to be negative, as in most materials used 
in microfabrication (glass, silicones, silica, etc.). We shall 
see that the surface promotes OLC only if its charge has 
the same sign as the inactive species (anions). 

Surface Conduction. (SC) — As the channel thick- 
ness is decreased, convection eventually becomes negli- 
gible compared to diffusion (small Peclet number), so 
we begin by considering only "surface conduction" asso- 
ciated with the excess counterfoils (active cations) that 
screen the wall charge [TBI, |2"D] , neglecting streaming cur- 
rent (which we analyzed separately |18j). For long, nar- 
row channels (H <C L), and thin double layers (Ad <C H, 
where Ad is the Debye screening length), the Nernst- 
Planck equations can be homogenized (area-averaged) as 
follows [HI HUE], 




0.6 
c 

0.4 



12 16 20 





fb) 


\ \ \p=20 




\ \v=40 \ 
V=60\ \ \ 





0.2 0.4 v o.e 



surface transport (SC, EOF) 

—WWW- 



reservoir o- 



-O membrane 



(c) 



bulk transport (CP) 



FIG. 2: Ovcrlimiting current in the surface conduction 
regime: (a) Dimensionless current- voltage relation, Eq. ( |34[ |, 
for p s — —0.01,-0.1 and (b) mean concentration profile of 
inert anions, Eq. (33 1, for p s = —0.1. (c) The equivalent 
circuit: an ideal diode for bulk CP in parallel with a constant 
shunt resistance for SC (or EOF, as in Fig. 3). 
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where c+ and c_ are the mean concentrations of cations 
and anions, respectively, D is the ionic diffusivity (for 
simplicity, the same for both species), r/> is dimensionless 
potential scaled to the thermal voltage, kT/ze. Equation 
(31 1 relates the cation flux to the current density j in 
steady state for a perfect cation-selective surface. Equa- 



tion ( 32 ) is a mean electroneutrality condition including 
both ionic and fixed surface charges, where p s = 2o~ s /H 
is the negative volume density of fixed charge. 



The homogenized ID model can be solved analytically: 

= ln(c), c - p s ln(c) = 1 - jx, (3) 
j = 1 - e~* - p s V, (SC) (4) 



where p s 



2zeco 



J 



is dimensionless fixed-charge density, 

,,.. ilh . is the dimensionless current density scaled 
by the limiting current, corresponding to the case with 
neutral side walls; c = c/cq, x = x/L, and voltage 
V = zeV/kT = —(f)(1) is the potential drop across the 
electrolyte, relative to 0(0) = at the reservoir. The 
current-voltage relationship, Eq. (4), is shown in Fig. 
2(a) demonstrating nearly constant conductance in the 
OLC regime. For small surface charge and/or wide chan- 
nels, < — p s <C 1, the anion concentration profile (33) is 
linear in the quasi-neutral bulk region, 5 ~ 1- 



cays exponentially, c 



,{jx-l)/p B 



-jx, and dc- 
in the depleted zone, 



j < x < 1, that forms above the limiting current, as 
shown in Fig. 2(b). In the opposite "membrane limit", 



— p s ^> 1, the concentration is only weakly perturbed by 
the current, c ~ 1 + (jx — l)/p s - 



As shown in Fig. 2(c), the current-voltage relation (34) 
can be interpreted as an ideal diode, j\> = 1 — e~ v , for 
bulk CP in parallel with a shunt resistance, j s = —p s V, 
for SC. The electric field, which is nearly uniform across 
the microchannel, acts on two different types of ions: (i) 
the inert anions plus an equal number of cations ("bulk 
conductivity" ) , and (ii) the remaining cations that screen 
the surface charge ("surface conductivity"). The re- 
sponse is equivalent to that of an ideal Schottky diode 
(metal/semiconductor junction) [T!5] with the surface 
charge playing the role of a n-type dopant in the semi- 
conductor, which provides a residual conductivity under 
reverse bias. As the microchannel thickness is increased, 
convection due to EOF eventually dominates SC, and the 
semiconductor analogy breaks down. 

Electro-osmotic Flow (EOF). — For our "dead-end 
channel" with an impermeable membrane, pressure- 
driven back flow opposes EOF and results in two counter- 
rotating vortices. Concentrated solution flows to the 
membrane along the side walls, and depleted solution 
returns in the center (Fig. 1(c)). 

As a first approximation, we assume thin double layers 
and a neutral bulk solution (c = c_ = c+) described by 
the steady 2D Nernst-Planck-Stokes equations, 

u-Vc = -V-F ± , F± = -D (VcicV^) (5) 
Vp = ?7V 2 u, V • u = (6) 

where u is the fluid velocity, 77 the viscosity, and p the 
pressure with boundary conditions of no slip, u = 0, and 
vanishing normal co-ion flux and specify uniform electro- 



FIG. 3: 2D simulations of the EOF mechanism for OLC (A = 0.01, H = 0.2, a a = —0.01): (a) Bulk salt concentration c 
(contours) and velocity vectors u (arrows); (b) Area-averaged concentration profiles (continuous), showing a nearly constant 
depleted-zone concentration Cd ~ l c - Inset: vortex position l c vs. voltage V; (c) Current-voltage relation from simulations 
(continuous) and the ID model (|8| (dashed). Inset: Area-averaged concentration profile in the 2D simulations (continuous) 
and in ID model (dashed) for V = 50 . 



chemical potential of counter-ions at the membrane cor- 
responding to potentio-static operating conditions; uni- 
form concentration c(0, y) — cq and zero pressure at the 
reservoir; and zero normal flux (neglecting SC) and EOF 
slip at y = ±H/2 given by the Helmholtz-Smoluchowski 
formula (neglecting diffusio-osmosis [3]), 
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where C( c ) is the ze ta potential in the Debye-Hiickel 
approximation and e, the permittivity. The main out- 
put parameter in model (|35|-([7|, the dimensionless cur- 
rent, j, depends on the dimensionless channel thickness, 
H = j^, reservoir Debye length, Ao = Ap ^ c °) , surface 
charge, a s — p s H, and voltage, V = 

Numerical solutions [18] of ([35])- (fTj) are shown in Fig. 
3. An interesting observation about CP during OLC fol- 
lows from averaging the 2D concentration profile (Fig. 
3a) over the channel cross section (Fig. 3b) . The intense 
vortical circulation near the dead end forms an extended 
depletion zone with a nearly constant area-averaged con- 
centration, Cd = CdCo, as the case of electro-osmotic in- 
stability (EOI) without geometrical confinement [3]. For 
moderate aspect ratios, < 10, the dimensionless dis- 
tance of the vortex center from the membrane, l c = l/L, 
is nearly independent of voltage in the OLC regime (inset 
of Fig. 3b) and approximately equal to the dimensionless 
concentration, dd ~ Re- 
using this observation, we can derive a ID model for 
the EOF regime [T5] (inset of Fig. 3c) which yields an 
approximate OLC current-voltage relation, 



which compares well with our 2D simulations (Fig. 3c) 
for H > 0.1 [IS]. The linear dependence of (18j again 
implies a constant shunt resistance (Fig. 2c) for the EOF 
mechanism, j s = CdV ', as in (34) for the SC mechanism, 



but with a different over-limiting conductance, Cd rather 
than — p s . 

To understand the over-limiting conductance, we per- 
form a scaling analysis of the 2D model [18J . Balancing 
axial convection over the length scale l c with transverse 
diffusion over the scale H in (35), the vortex position 
scales as, l c ~ u E q F H 2 /D. The characteristic EOF ve- 
locity along the walls scales as, ueof ~ &sju re f /c 3 / 2 \o, 
where we use ([7| and define the reference (electroviscous) 
velocity scale, u re f — e(kT/e) 2 /r)L. The numerical ob- 
servation, Cd ~ l c , then implies the following scaling re- 
lation near limiting current (j s» 1): 



Cd ~ l c 



0.22a s 2 / 5 ff 4 / 5 A - 2/5 , 



(9) 



j w 1 - cj + Tc d In Cd + CdV (EOF) 



(8) 



which shows how the over-limiting conductance due to 
EOF varies with surface charge, channel thickness, and 
reservoir concentration (via Ao). Our 2D simulations 
(Fig. 4, inset) roughly confirm the scaling of Cd ~ l c 
with H and provide the prefactor in For very large 
currents (j 3> 1), the conductance is predicted to in- 
crease [IB] (Fig. 3c), but a detailed analysis also con- 
sidering high aspect ratios (H^ 1 ^> 10) is left for future 
work. 

Unified scaling analysis. — The relative importance 
of the SC and EOF mechanisms for OLC is controlled 
by the dimensionless ratio of their conductances in the 
weakly over- limiting regime, £ = \p s \/cd- Using (JoJ) and 
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FIG. 4: The regions of dominance of the SC, EOF, and EOI 
mechanisms for varying channel aspect ratio, H~ and surface 
charge a s for Ao = 10~ 5 . The critical value of H is given by 
£ = 1 in jioj) for the SC/EOF transition and by || = in 
2D simulations for the EOF/EOI transition. Inset: Simulated 
vortex center location l c versus H for Ao = 0.01, a s — —0.1 
(solid) compared to the scaling law, f c = 0.38JT 1 / 5 in J9J. 



restoring dimensional variables, we obtain the scaling, 
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where S = 1 marks the transition between mechanisms 



(Fig. 4). Equation (10) shows that SC dominates for 
thin channels, low ion concentration and/or large surface 
charge (S 3> 1), while EOF dominates above a critical 



channel thickness and/or salt concentration (S <C 1). 

As the channel thickness is increased further, sur- 
face effects eventually become negligible, and OLC re- 
sults from bulk EOI at the membrane, as in unconfined 
systems [3J 2]. The EOF vortex intensity, scaling as 
l c /H ~ iJ^ 1 / 5 from (lij, decreases with H, and thus, 
above a critical channel thickness, a significant portion 
of the depleted solution near the membrane is unaffected 
by EOF and instead undergoes EOI to maintain OLC. 
EOI also leads to a constant over-limiting conductance 
(beyond a threshold voltage) [1 H], so the EOF/EOI 
transition can be defined via 4% = in our 2D sim- 
ulations, above which the EOF supply of solute to the 
membrane fades and conditions develop for EOI (Figs. 
4 and S4 [18]). In water, for L = 1mm, ( = 50m V 
and Co = ImM, the predicted mechanism for OLC tran- 
sitions from SC to EOF at H = 8/im (consistent with 
observed vortex supression [IT]) and from EOF to EOI 
at H = 0.4mm, so that different microfluidic experiments 
can be dominated by SC [Hj, EOF [JO], or EOI [3J. 

Conclusion. — In a microchannel, charged side walls 
can sustain OLC to an ion-selective end. With increas- 
ing channel thickness, the dominant mechanism switches 
from SC to EOF to EOI. These nonlinear phenomena 
may find applications in water desalination, biological 
separations, or electrochemical energy storage. 
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Supplementary Information 



1. EFFECT OF STREAMING CURRENT ON 
SURFACE CONDUCTION 

Equations (1) and (2) in the main text present a model 
for the area-averaged concentration in a channel. In this 
model axial transport of ions is considered via diffusion 
and electromigration including migration of excess ions 
shielding the wall charge, but the effects of flow is ignored 
(streaming current). We here present an estimate of this 
effect in the simplest form and show that consistent with 
our assumption, this effect is negligible for sufficiently 
thin dead end channels. 

It should be noted that the relative importance of 
streaming current in surface conduction in the context 
of thin double layer is well established in the literature. 
For example, Deryagin and Dukhin 19 used the Gouy- 
Chapman EDL model and derived analytical expressions 
for streaming current carried by electroosmotic flow. Fol- 
lowing a similar analysis one could show that the ratio 
of the streaming current to current carried by electro mi- 
gration of excess counterions shielding the wall charge 
is: 



Jilo 



7cond. 



2Pe tanh 



4kTj ' 



(11) 



where Pc = e(KT/ez) 2 jr\D is the material Peclet num- 
ber and property of the solution. This analysis however, 
is not applicable to the case of channels considered in 
our studies; when channels themselves are thin, the thin 
EDL assumption will be violated due to overlap EDL ef- 
fects. In addition, the pressure driven flow, which would 
cancel the net electroosmotic flow for dead-end channels, 
will contribute to current. 

To remedy these issues we carried out a separate anal- 
ysis specific to the case of dead-end thin channels. One 
can consider a channel with high aspect ratio (L 3> H); 
therefore the competing flow components, i.e. pressure 
driven and electroosmotic flow, can be assumed to be 
locally parallel. Furthermore, the EDLs (either over- 
lapped or non-overlapped) can be assumed to be in local 
equilibrium 12 . Therefore, the charge distribution can be 
estimated from solution of the Poisson-Boltzmann equa- 
tion with top and bottom wall condition as known surface 
charge, a s . These solutions determine the distribution 
of the potential and charge in the wall-normal direction. 
Having such solutions one can determine the velocity pro- 
file associated with the electro-osmotic flow: 



eE 
V 



Avail 



(12) 



where E is the axial electric field, and cj>(y) — </> W aii is 
the local potential difference in the wall-normal direction. 
The pressure driven flow is assumed to have a parabolic 
profile with a vanishing net flow. In such a scenario, near 



the walls the electro-osmotic flow dominates the pres- 
sure driven flow and advects the charged double layers 
forward; in the middle section, the pressure driven flow 
wins, but the available charge is much less than that near 
the wall. Therefore, the flow would contribute a net cur- 
rent to the system. This net current is proportional to 
the flow, which itself is proportional to the applied field, 
and thus has an Ohmic effect. 

This problem can be solved in the general form. 
The two dimensionless parameters associated with the 
Poisson-Boltzmann system are: Xjj/H, Debye length 
over channel thickness, and p s , the dimensionless wall 
charge. The current contribution of the flow, becomes 
more important in the limit of highly depleted channel, 
i.e. Xd/H — > oo, and p s — > oo. In this limit, the mean- 
ingful controlling parameter is 



aHez p s H 2 



2eKT 



4X1 



(13) 



which is finite and independent of the bulk concentration. 
Figure [5] shows the net streaming current as a function 
of a* from our analysis and contrasts it to the prediction 
of thin EDL model (Eq. 11), which is clearly shown to 
overestimate the effect of streaming current. 

For a one micron thin channel with surface charge 
<7=lmC, er* will be about 30; assuming an aqueous solu- 
tion at room temperature (Pe~0.5) the current contribu- 
tion by the flow is estimated to be about 15% of current 
by the electromigration of the surface charge. In this 
case ignoring the streaming current is marginally valid, 
and the approximation significantly improves as one uses 
thinner channels, however, for larger channels (or larger 
surface charge) the streaming current should be consid- 
ered in the model (in simplest form, as a prefactor en- 
hancing the surface conduction effect). 



2. EOF MECHANISM: SCALING ANALYSIS 

Here we derive in more detail the scaling relation 
Eq.(9) for the over-limiting conductance due to electro- 
osmotic flow on the side- walls of the microchannel, which 
leads to a pressure-driven backflow from the imperme- 
able membrane end of the microchannel. The resulting 
vortex pair drives the fast convection of salty solution 
to the membrane along the walls and the return of de- 
pleted solution back into the microchannel (Fig. 2c). It 
is tempting to analyze convection-diffusion in this flow by 
area-averaging using the theory of Taylor- Aris dispersion, 
but we shall see that convective mixing in the transverse 
direction is not fully developed in the vortex region and 
leads to different scalings with system parameters. 

Electro-osmotic flow near the channel wall scales as 



(14) 



where the zeta potential, C, is related to the surface 
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FIG. 5: Streaming current versus nondimensional wall charge. 
Solid line represent the result of our analysis for thin chan- 
nels including effect of the pressure driven flow. Dashed line 
represents the thin-EDL approximation based on Eq. |11| for 
X/H = 0.5. Current is normalized by the current associ- 
ated with electromigration of the surface shielding charge, 
ezv\a\E (per channel width). The plotted result is valid for 
all Xd/H > 0.5 within plotting accuracy. 



charge density using the thin diffuse-layer capacitance 
(C = e/Xd), and the tangential potential gradient, (j) x , 
is related to the current density, j: 



c 



ct s Xd 

e 



J_ 

Kb 



(15) 



where Xd(c) is the Debye length and Kb(c) is the bulk 
conductivity, each depending on the local bulk salt con- 
centration c. 

Combining these equations, we obtain the velocity scal- 
ing 



(16) 



Near the dead end, the flow recirculates and produces 
a vortex which shears the fluid over a transverse scale 
H (the channel thickness) and longitudinal scale l c (the 
vortex center). The concentration is strongly depleted 
by the selective surface (e.g. membrane) near the end of 
the channel, where the electric field driving the flow is at 
maximum. Near this end the flow turns and, upon mov- 
ing away from the end, eventually develops into a plane 
parallel slip/pressure driven flow for which the vortex 
center location is similar to an entrance length. This re- 
sult is reminiscent of the convection-diffusion problem for 
the evaluation of the entrance length of a cold fluid en- 
tering a hot parallel-plate pipe. In our terms, the length 
of the diffusion scale is the channel thickness H, while 
the longitudinal convection scale is the vortex center lo- 
cation (the entrance length) l c . The convection-diffusion 
equation, u ■ Vc = £>V 2 c, yields the scaling 



D 
H 2 



(17) 



which implies 



uH 2 cr s X D jH 2 



D 



K b ?]D 



(18) 



Next, we introduce dimensionless quantities, denoted 
by tilde accents. Let us employ the channel length L for 
the length scale, kT/ze for the potential, and thus the 
electro- viscous scale, 



Iref 



sjkT/zef 



(19) 



for the velocity scale. Let us also scale concentration 
to the bulk reservoir concentration Co at the open end 
of the channel. Assume an z : z electrolyte with equal 
diffusivities and an ideally cation selective surface. This 
yields the following expression for conductivity and De- 
bye length, 



K b 



eD 2 _ ekT 



2(ze) 2 c 



(20) 



The dimensionless values for current, height and surface 
charge, are defined as 



J 



JL 



2zeDcn 



L ' 



zecoL 



(21) 



This nondimensionalization yields two crucial parameters 
the "material Peclet number", Pe, (as 1/2 for aqueous 
solutions at room temperature) and dimensionless Debye 
length, Aq, at the entrance to the channel, defined as 



Pe 



u re fL e{kT / ze) 



D 



\ 2 

A 



ekT 



2(ze) 2 c L 2 



(22) 



With these definitions, the scaling of the vortex center 
takes the dimensionless form 



Pe H 2 a s j 
A 5 3 / 2 



(23) 



The only unknown in the scaling ( 23 ) is the dimension- 
less concentration in the mixed depleted region. Using 
the numerically determined relationship 



lc 



(24) 



which is valid in the weakly overlimiting regime (Section 
2 below), we obtain 



PeH 2 a s 
A 



2/5 



(25) 



where the coefficient of proportionality is approximately 
0.3 from simulations. As this equation is applicable only 
in the weakly overlimiting regime (j ~ 1), the current 
density does not affect the scaling of the overlimiting 
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conductance. If we were to keep the current from our 



(a) 



(b) 



derivation and still use ( 24 ) , then we would obtain 

2/5 



Cd 



Pe H 2 d s j 
Ao 



(26) 



for very large currents, j> 1, in place of (25) 



To compare the EOF vortex mechanism for over- 
limiting current (dominant for wide channels) to the 
mechanism of surface conduction (dominant in narrow 
channels) in the weakly overlimiting regime, we consider 
the dimensionless parameter E, 



fh 



~c d H 



(27) 



which is the ratio of the two different possible depth- 
averaged conductivities in the depleted region, where 



Pa 



zecoH 



1l 

H 



(28) 



is the scaling of the limiting conductivity due to surface 
conduction. Using Eqs. ( 24 )- ( 25 1 and the simulation re- 
sults, we obtain 



Pe 2 H 9 



1/5 



(29) 



which depends on the channel dimensions, surface charge 
density, and reservoir salt concentration as 



H 



-9/5 t-4/5 _3/5 -4/5 



For sufficiently narrow and long channels, thick double 
layers, large surface charge density, and/or small reser- 
voir concentration, £ 3> 1, surface conduction is the 
dominant mechanism for over- limiting current. In the 
opposite limit, S< 1, the dominant mechanism is con- 
vective mixing by confined electro-osmotic flow. 



2D SIMULATIONS 



Here we study in detail the EOF mechanism of the 
OLC in the dead-end microchannel via 2D numerical 
simulations in the limit of thin double layers. The rel- 
evant set of steady-state dimensionless Nernst-Planck- 
Stokes equations for the electroneutral bulk, neglecting 
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FIG. 6: (a) Electric potential, <p, level lines; (b) concentration, 
c, level lines and (c) EOF streamlines computed for N = 10, 
A = 0.01, cr s = 0.01, V = 50 and H = 0.2. 



SC, reads 



Pcu • V2 = V 2 a, V(cV<£) = 0, V 2 u - Vp = 0; (31) 

V = 0, H : ^ = ^= 0,u-n = 0, u-r = cfe (32) 
oy oy ox 

3c 3(b 

x = l; — -cV- =0,u = 0, ln5 + (^ = lniV; (33) 
ox ox 

d 

x = 0: c = 1, <p = V, u • t = 0, — (u • n) = 0. (34) 

Ox 



(30) Here £ in the Helmholtz-Smoluchowsky slip condition 



(32 1 is the dimensionless ^-potential, which we assume is 
related to the surface charge of the nonconducting side 
wall through the linear (low- voltage) Debye Huckcl rela- 
tion, 



c 



(35) 



neglecting any effects of the compact Stern layer. 
Conditions (33) at the perfect cation-selective mem- 



brane/solution interface are those of vanishing anions' 
flux and continuity of the electrochemical potential of 
cations across the interface, along with the non-slip con- 
dition; N is the dimensionless fixed charge density in the 
membrane and V is once more the potential drop across 
the diffusion layer. 

Below we present the results of numerical solution of 
the problem (31-34). Thus, in Fig. 2 we present the 



concentration and electric potential level lines and EOF 
streamlines obtained for highly OLC regime. In Fig. 3 we 
show the results of numerical verification of scaling re- 
lations (Eq.(9), main text) derived in the main part of 
this Letter. Thus, in Fig. 3a we present the numerically 
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FIG. 7: Simulation of EOF mechanism for Ao = 0.01: (a) The coefficient of the proportionality Iq versus channel depth H 
for j — 1.5 and the following values of the surface charge density: (1) - a s — 0.02, (2) - a s = 0.025, (3) - a s = 0.03, (4) - 
o s = 0.035, (5) - a s = 0.04; (b) Concentration of the depleted bulk Cd versus the electric current density j for H = 0.2 and 
o s = 0.04, compared with j' 1 ^ 3 and j 2//5 asymptotic power laws (dashed lines); (c) the ratio Cd/l c versus surface charge density 
cr s computed for j = 1.5 and H = 0.2. 



computed coefficient of proportionality Iq 



(36) 



versus channel depth H for a sequence of values of sur- 
face charge density d s . In Fig. 3b we show the weak de- 
pendence of the depleted bulk concentration Cd on the 
electric current density in the weakly overlimiting cur- 
rent regime, which transitions to power-law dependence 
in the strongly over-limiting regime. Although the simu- 
lations may be consistent with an asymptotic power law, 
Cd ~ 7 2 / 5 , as predicted by our simple scaling analysis, 
Eq. (26 1, a somewhat weaker dependence, q ~ j 1 ^ 3 



provides a better fit over the range of moderately large 
currents in the simulations (up to j = 6). In Fig. 3c we 
show that the increasing rate of the EOF due to an in- 
crease of surface charge density yields convergence of the 
dimensionless concentration of the depleted bulk Cd to 
the dimensionless location of the vortex center l c for the 
weakly overlimiting current regime. 

Finally, to illustrate the transition to EOI mediated 
OLC with the decrease of the channel aspect ratio in 
Fig. 4a we present plot of the dimensionless concentra- 
tion of the depleted bulk Cd versus H which shows a sharp 
decrease of Cd for H above some threshold due to the in- 
creasing insensitivity of the middle part of the channel to 
the wall-induced EOF. In Fig. 4b, c we illustrate this by 
presenting the EOF streamline maps for H = 0.2 below 
and H = 1.6, above this threshold. 



4. APPROXIMATE ID MODEL 

Motivated by our 2D simulations, we can develop a 
simple ID phenomenological model for the EOF driven 
OLC in which l c is the sole fitting parameter. In this 
model the dimensionless salt concentration in the dc- 




FIG. 8: (a) The dimensionless concentration of the depleted 
bulk c d versus H for a s = 10" 3 and V = 10; (b) EOF stream- 
lines for o s = 10~ 3 , V = 10 and H = 0.2; (c) EOF streamlines 
for d a = 10" 3 , V = 10 and H = 1.6. 
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pleted region is fixed at Cd = l c , whereas in the rest of 
the diffusion layer it varies linearly from l c to unity at 
the open end of the channel and from l c to l\ at the dead 
end of the channel: 



C let 



1 - L 



< x < 1 



L - li 



J J 
i - L 



c = 1 — jx, < x < 



c = P c +j(l-x), 1 



J 

I - ll 



<x<l: 



The resulting current /voltage relation 



lc 



C 1 ~ 5 



(37) 
(38) 
(39) 



c-^ = -j, 0<x< 1,0(1) -<^(0) = -K (40) 



(41) 



is plotted in Fig. 3c (main text) assuming constant Cd ~ 
l c . Comparison with the 2D simulations allows us to 
roughly confirm the scaling relation ( 25 ) and fit the pref- 
actor. 



We can also begin to understand the appearance of 
a weak nonlinearity in the 2D simulation results (Fig . 
3c, main text), where the over-limiting conductance due 
to EOF grows slowly with voltage. Taking into account 
the weak variation in Cd for j from the scaling analysis, 
we would predict a weak nonlinearity, j ~ t/ 5 / 3 for very 
large currents, j>l. Alternatively, assuming Cd ~ j 1 ^ 3 
from Fig. 2b would suggest j ~ V 4 ^ 3 , which may also be 
consistent with the simulations. A detailed analysis of 
the strongly over-limiting regime is left for future work. 



